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Abstract 

We  discuss  the  interaction  of  a charged  particle  with  the  surface  excitations 
of  a semi-infinite  medium.  We  consider  only  the  case  where  the  particle  is 
outside  the  material  and  trapped  near  the  surface  by  its  image  potential.  By 
the  use  of  a cononical  transformation  which  exhibits  the  bound  states  due  to 
the  image  charge  in  the  non-perturbed  transformed  Hamiltonian,  we  determine 
the  binding  energy  and  the  effective  mass  of  a charged  particle  for  motion 
parallel  to  the  surface.  For  the  highest  perpendicular  sub-bands  the  correc- 
tions to  the  energy  and  mass  of  the  particle  are  negligible.  We  have  applied 
the  method  to  determine  the  binding  energy  of  the  ground-state  of  an  electron 
trapped  near  the  LiF  surface.  We  find  the  value  E(0)  “ -0.25  eV. 
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I.  INTRODUCTION 


Both  theory  and  experiment  have  been  extensively  developed  by  treating 
the  interaction  of  charged  particles  with  collective  excitations  on  surface 
media;  for  example,  surface  plasmons  in  metals  and  surface  optical  phonons  in 
dielectric  crystals.^  This  interaction  with  surface  modes  can  occur  when  the 
particle  is  either  inside  or  outside  the  surface.  We  consider  here  the  case 
where  the  charged  particle  lies  outside  the  material,  interacting  then  only 
with  surface  modes.  When  the  particle  is  inside  the  medium  it  is  also  necessary 
to  treat  the  interactions  with  the  bulk  modes. 

It  is  well  known  that  one  charged  particle  outside  the  material  induces 
charges  of  opposite  sign  on  the  surface  and  is  attracted  to  it  by  the  image 
potential.  It  can  form  a bound  state  in  which  it  is  trapped  near  the  surface 
of  the  material  by  the  attractive  image  charge.  It  is  also  well  known  that 
the  image  potential  is  due  to  the  surface  modes  of  oscillations.  These  modes, 
when  polarized,  create  electric  fields  which  act  back  upon  the  electric  charge. 

This  is  the  source  of  the  image  potential. 
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Cole  and  Cohen,  and  Cole  have  showed  that  surface  states  induced  by  the 
electrostatic  image  potential  at  insulators  such  as  He,  Ne,  H2  and  D2,  both 
solid  and  liquid,  have  too  small  binding  energies.  In  these  weakly  bound  states, 
the  electron  is  localized  well  away  from  the  surface. 

In  this  work  we  investigate  the  possibility  that  a charged  particle  may 
be  trapped  near  the  surface  of  a semi-infinite  material  by  the  charge-surface 
excitation  interaction.  The  material  may,  in  principle,  be  either  solid  or 
liquid.  Our  major  purpose  is  to  determine  the  binding  energy  and  the  effective 
mass  of  the  ground  state  of  a charge  particle  placed  above  the  surface  of  a 
material . 

We  shall  assume  that  the  difference  between  the  vacuum  level  and  the  bottom 
of  the  conduction  band  is  much  larger  than  the  binding  energy.  In  this  situa- 
tion we  may  simplify  the  problem  by  considering  that  the  charge  cannot  penetrate 

the  material.  This  approximation  is  discussed  in  detail  both  in  reference  (3) 

4 

and  in  a recent  review  paper  by  Crandall  for  the  case  of  electrons  trapped 

4 

near  the  surface  of  He  . It  is  probably  also  applicable  for  LiF  where  the 
potential  barrier  is  much  larger  than  the  binding  energy.^  In  these  cases,  the  elec 
tron  wavefunction  is  concentrated  in  a region  far  from  the  surface  and  the  energy 
required  for  the  electron  to  penetrate  the  material  is  too  high.  We  also  assume 
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that  the  interaction  between  the  charged  particle  and  the  surface  excitations 
is  very  weak,  which  means  that  the  particle  will  be  loosely  bound  to  the  surface. 

We  proceed  in  a similar  fashion  to  the  canonical  transformation  technique 
employed  by  Freitas  and  Lobo.^  As  it  was  pointed  out  in  reference  (6),  the 
canonical  transformation  chosen  has  the  merit  of  exhibiting  the  bound  states 
due  to  the  image  potential  in  the  unperturbed  part  of  the  transformed  Hamiltonian. 

In  Section  II  we  describe  the  model  to  be  employed  and  we  discuss  the 
interaction  between  a charged  particle  and  longwavelength  excitation  modes  in 
a semi-infinite  material.  We  show  in  the  framework  of  classical  electrodynamics 
that  a charged  particle  placed  outside  the  material  can  excite  only  surface 
modes.  We  employ  a canonical  tmasformation  technique  for  the  whole  Hamiltonian 
of  the  system  and,  in  the  limit  of  weak-coupling  constant,  part  of  the  kinetic 
energy  of  the  particle  is  treated  as  a perturbation. 

After  this  examination  of  the  interaction  of  the  charged  particle  with 
the  surface  modes,  we  present  in  Section  III  the  results  for  the  energy  shift 
which  arises  in  the  second-order  of  perturbation  theory.  We  calculate  the 
binding  energy  of  the  ground-state  of  a particle  trapped  near  the  surface  as 
a function  of  , the  coupling -constant,  and  also  the  effective  mass  of  the 
particle,  for  motion  parallel  to  the  surface. 


II.  MODEL  HAMILTONIAN 


Consider  a charged  particle  placed  a distance  z>0  from  the  surface  of  a 
semi-infinite  material.  We  take  the  material  to  fill  the  half-space  zKO, 

We  assume  that  the  deielactric  constant  s (^)  is  frequency-dependent  only.  A 
more  general  description  would  have  the  dielectric  constant  frequency  and  wave- 
vector  dependent.  The  approximation  of  neglecting  the  wave-vector  dependence 
gives  equations  which  are  much  easier  to  solve  although  it  has  no  justification 
in  all  physical  situations. 

For  the  sake  of  clarity  and  completeness,  we  shall  show,  in  a standard 

procedure,  that  in  the  longwavelength  limit  an  external  charged  particle  does 
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not  interact  with  bulk  modes.  ’ The  reason  is  that  the  fields  associated  with 
bulk  modes  are  confined  to  the  material  only.  The  interaction  between  an  elec- 
tric charge  e and  a given  collective  mode  is  given  by  e^(r),  where  fl(r)  is 
the  electrostatic  potential  associated  with  this  oscillation.  Since  there  is 
translational  invariance  along  a direction  parallel  to  the  surface,  the  type  of 
solutions  we  seek  oarrespond  to  waves  for  which 


(r,t) 


'S(Z) 


-idUt 

e 


(2.1) 


where  R is  the  projection  of  the  position  vector  r onto  a plane  parallel  to 
the  surface  and  Q is  a two  dimensional  wave-vector. 

Assuming  that  there  is  no  charge  density  or  currents  in  this  system,  we 
have  from  Maxwell's  equations 


e(uj) 


z < 0 


and 


(2.2) 


dz 


z > 0 


(2.3) 


The  solutions  of  Eq.  (2.2)  are  either  s (a))  =0,  or  a linear  combination 


of  the  two  Independent  solutions,  e*^^  and  e 


-Qz 


The  solution  s (o)  “ 0 gives 


the  frequency  of  the  bulk  excitations  and  the  second  possibility  describes  the 
surface  modes.  Outside  the  material  v(z)  satisfies  Eq.  (2.3)  and  is  given  by 


5(z)  = e 


(2.4) 


and  for  the  bulk  inodes  solution,  -(z)  assumes  the  form 


cp(z)  = V sin  (q  2)  + -p:,  cos  (q  z), 


(2.5) 


where 


^1’  ^2  ^3  constants  to  be  determined. 


The  boundary  conditions  at  the  surface,  z = 0 , say  that  the  normal 
components  of  D and  the  tangentical  components  of  E must  be  continuous. 

Since  € (u))  = 0,  we  have  = 0 inside  the  material.  Then,  just  outside 

D = 0 , which  requires  from  Eq.  (2.4),  t)-  s 0 . So,  the  electrostatic  poten- 
tial  ^(r)  is  identically  zero  outside  the  medium  which  means  that  the  parallel 
electric  field  is  zero  just  outside.  The  continuity  of  the  parallel  components 
of  E at  z = 0 requires  from  Eq.  (2.5),  ^ 0 .• 

We  conclude  that  the  electrostatic  fields  associated  with  bulk  oscilla- 
tions are  confined  inside  the  material  where  the  charge  oscillations  take  place. 
These  fields  are  derived  from  an  electrostatic  potential  d(r,t)  , given  by 
the  following  expression 


- iQ.R-i(JJt  . , . 

I(r,t)  - r-e  ^ sin  (q  z) 

i.  z 


(2.6) 


By  examining  the  solutions  of  Eq.  (2.2)  which  describe  the  surface  modes 
and  by  matching  at  the  surface,  the  fields  inside  and  outside  of  the  material, 
we  obtain  the  condition 

e (ID)  = - 1 

which  determines  the  frequency  of  the  surface  oscillation  modes. 

Then  the  electrostatic  fields  associated  with  these  modes  are  extended  outside 
the  material,  an  important  distinction  with  bulk  modes.  In  other  words,  a charged 
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particle  outside  the  system  can  excite  only  surface  modes. 

In  this  work  we  shall  discuss  only  the  interaction  of  the  low-energy 
charged  particle  with  surface  modes.  By  assuming,  for  simplicity,  chat  the 
charge  cannot  enter  the  material,  we  may  write  the  Hamiltonian  of  the  system  in 
three  parts  : 


« ' «C  ^ «S  ^ »int  « 


(2.7) 


where  is  the  charged  particle  kinetic  energy,  Hg  describes  the  surface 

modes  in  the  absence  of  the  particle  and  is  the  interaction  between  the 

charged  particle  and  the  surface  modes.  The  charged  particle  part  has  the  form 


= - Qr~}/2m)  7^  , 


(2.8) 


where  m is  the  mass  of  the  particle.  The  surface  modes  part  of  the  Hamiltonian 


can  be  written  as 

Q 

where  the  surface-wave  frequency  is  given  by  the  relation  ®('^g)  ~ 
a^  and  a^  are  the  creation  and  annihilation  operators  respectively  for  the 
surface  excitations,  obeying  the  usual  commutation  relations, 


(2.9) 


[®Q  ’ I ~ 'QQ* 


(2,10a) 


["q  ’ "q’1  • 4J  = ° 


(2.10b) 


The  interaction  between  the  particle  and  the  surface  modes  is  given  by  the 
, 1 

expression  ; 


H,  ^ = 2-1 
int 


(£r)l 


(2.11) 


where  A is  Che  area  of  the  surface,  z is  the  distance  of  the  particle  from 
the  surface  and  the  parameter  5 has  different  value  for  each  type  of  excita- 
tion. For  example,  it  will  be 


3 = uUg  /2tt 


3.!§_  /Vi  . !£i) 


plasmons 
phonons , 


(2,12a) 

(2.12b) 


with  s„  , the  static  dielectric  constant  of  the  material  and  s , the  optical 
S 


dielectric  constant. 

In  analogy  with  the  bulk  excitations,  we  define  the  coupling -cons tant 
as 


(*) 


(2.13) 


In  a manner  similar  to  Freitas  and  Lobo,  we  introduce  a canonical  trans- 
formation S which  removes  the  coordinate  R and  gives  a convenient  zeroth 
order  Hamiltonian  with  bound  states  for  the  external  charged  particle. 

We  define  the  unitary  operator 


S = exp 


*5  - f 


5 *0 


(2.14) 


with 


-iQ.R-Qz 


f-  = - (2ni)  (e/ha>g^  (hp/2AuUg) 


(2.15) 


The  new  Hamiltonian  H will  be  generated  by  application  of  the  unitary 
transformation  S in  the  standard  manner: 


y = S"^  H S 


(2.16) 


The  explicit  form  of  ^ after  expanding  up  to  second  order  over  the  coupling 
constant  has  the  following  expression; 
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y = - ^ 3-  3- 

2m  ^ 


+ 

S ”Q  "Q 


-A 


I 2rr\'?^ 


mA 


-Qz 


where 


• X 

2m 


r r o' 

n 

) 

[u. 

+ ^!  U,  U,  T 

+ . • • . • 1 

) 

u =E 


L ('q  *q  - £q  >5  ). 


and 


is  a commutator. 


(2.17) 


(2.18) 


[•] 

By  performing  the  Q summation  on  the  third  term  of  the  right-hand  side 

2 2 2 

of  Eq.  (2.17)  and  by  combining  this  result  with  the  3 /dz  term  of  .we 
get 


1 

.2  ,2  \ 

2 + 2 1 

- K 

n ? e^  \ 

2m  \ 

. 2 2 I 

ox  oy  / 

^ 2m 

""s'  " ' 

Q 

2m 

^ ! 

Q 

f 

1 

1 

* _2 

.....i 

1 

1 

. (2.19) 
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In  the  weak-coupling  limit,  first  order  in  OTg  , the  terms  which  come  from 
the  commutator  in  Eq . (2.19)  are  treated  as  perturbations  to  the  new  unperturbed 
Hamiltonian,  . 


We  write  the  unperturbed  state  vector  corresponding  to  the  total  momentum 
K parallel  to  the  surface  in  the  form 


1:>  = X_  (z)  ' V > , ( 

2-T  " q 'I 

where  are  the  normalized  s -wave function  of  the  one-dimens ional  hydrogen 

like  Hamiltonian, 


/ N - 2z 

~ 3/2  2 

a n 


(„!/  I “ / "I"*/ 


2 2 2 1 . 
where  a = (2iVg  .1  /rime  3)  is  the  effective  Bohr  radius,  is  a modified 

Laguerre  polinomial  and  qKq  ^ wavefunction  of  the  excitation. 

In  the  next  section  we  attempt  to  determine  the  two  quantities  o';  interest 

here,  namely,  the  binding  energy  E(0)  , and  the  effective  mass  m*  of  the 

ground-state  of  a charged  particle  for  motion  parallel  to  the  surface 


III.  BINDING  ENERGY  AND  EE’FECTIVE  MASS 


It  is  a straightforward  matter  to  calculate  the  energy  of  a charged 
particle  for  the  wavefunction  described  in  Eq.  (2.21).  Strictly  speaking, 
we  will  determine  the  energy  within  the  zero  plus  one  excitation  mode  space 
and  in  the  limit  of  weak-coupling  constant. 

In  this  approximation  the  eigensolutions  of  are 

u 

^n  • (3.1) 


where 


(3.2) 


is  the  energy  spectrum  of  the  s-hydrogen-like  atom.  Eq . (3.1)  describes  a set 
of  surface  sub-bands  with  band  index  n just  as  in  the  semiconductor  surface 
inversion  layers  problem.^ 

The  first  correction  to  the  energy,  particle  arises  in  the 

second  order  of  perturbation  theory.  Thus,  we  have 


n',K',p 


In  the  process  described  by  this  second-order  calculation,  the  intermediate 
states  correspond  to  a particle  in  state  n'  and  K'  and  to  one  surface  excita- 
tion with  wave-vector  p . 

It  is  very  difficult  to  get  values  of  the  correction  to  the  energy  and 
effective  mass  from  this  general  expression.  Therefore,  we  simplify  the  calcu- 
lations by  assuming  that  we  are  interested  only  in  the  case  where  the  particle 

remains  on  the  lowest  perpendicular  sub-band.  This  is  the  situation  discussed 
4 

by  Crandall  to  explain  the  long  life-time  of  an  electron  trapped  near  the 
4 

surface  of  He 


r.-.T’ 
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For  this  case  we  find 


(2)  - 2^  2 64 « 

El  ) 

" S 2Aa). 


-z 


/ .2,  2 .2.  2 \ 

I h K h P j 

\ 2m 2m  / 


p (2-^ap) 


h'"  - - 

— K.p  - 
m 


(3.4) 


Let  us  remind  the  reader  that  the  significant  contribution  to  the  sum  in  Eq. 

u 

(3.4)  comes  from  values  p < l/r,,  , r^.  “ (ll/2mu)  ) ^ a length  the  order  of 

the  radius  of  a bulk  excitation,  because  of  the  exponential  dependence  of  Eq . 
(2.11).  On  the  other  hand,  for  a 2r  , we  may  approximate  24^p  ^ 2 . By 
replacing  the  sum  by  integrations  according  to 


we  get  the  following  value  for  AE(0), 


AE(0)  = - 0.119  • 


This  result  fits  very  well  Evans  and  Mills'  theory  of  surface  polarons, 
provided  by  the  weak-coupling  and  no  penetration  of  the  variational  particle 
wavefunction  into  the  material.  As  can  be  seen  from  Eq . (3.5)  for  small  values 
of  Qfg  , the  particle  is  weakly  bound  to  the  surface  which  means  that  the  bind- 
ing energy  is  too  small  compared  to  the  surface  excitation  energy  hujg  . 

Strictly  speaking,  our  results  are  valid  only  for  smaller  than  1 but 

they  fit  quite  well  the  results  given  in  reference  (10)  for  the  binding  energy 

for  a as  large  as  4.  The  physical  picture  of  this  small  binding  energy  can 
o 

be  understood  in  the  following  way;  to  localize  the  particle  near  the  surface 
it  costs  a certain  amount  in  kinetic  energy.  But,  for  small  Og  , the  decrease 
in  energy  due  to  the  particle-surface  oscillation  modes  interaction  in  clothing 


the  particle  with  the  surface  excitation  is  nearly  cancelled  by  the  cost  in 
kinetic  energy. 

We  also  have  estimated  the  correction  to  the  ground-state  binding  energy 
which  arises  from  the  second  perpendicular  sub-band.  However,  we  found  that 
this  contribution  was  too  small. 

We  applied  this  veak-coupling  limit  method  to  estimate  the  binding  energy 
of  an  electron  trapped  near  an  ionic  crystal  surface.  For  example,  for  L^F 
the  coupling-constant  is  9fg=7.1  and  Eq.  (3,4)  gives  for  the  binding  energy 
of  the  ground -state  the  value 

E(0)  = -0.25  eV  , 


which  is  close  to  the  result  obtained  by  Sak.^^  Therefore,  we  feel  that  the 
weak-coupling  limit  is  not  realized  for  ionic  crystals  since  the  electron 
interacts  very  strongly  with  the  surface  oscillation  modes. 

We  now  report  the  effective  mass  results  of  the  particle  which  is  related 
to  the  coefficient  of  the  quadratic  term  in  the  expanded  energy,  as 


E(K)  = E(0)  + 


2m 


+ 


(3.6) 


Thus,  by  expanding  Eq.  (3.4)  in  terms  of  K , we  get  for  m the  expression 


* 1 
m - m (l+Tr  -T  O'.) 


16 


(3.7) 


As  can  be  compared,  the  correction  to  the  mass  of  the  particle  is  /l(s 

11  ^ 

whereas  it  is  for  the  weak-coupling  constant  theory  of  Sak. 

As  was  pointed  out  in  Sak's  paper  and  also  as  can  be  seen  from  Eq.  (3.4), 
the  renormalization  factor  of  the  mass  depends  on  the  functions  x^(z)  . It 
approaches  unity  as  the  energy  levels  of  the  perpendicular  sub-bands  become 
shallower.  Thus,  our  results  give  the  maximum  enhancement  within  the  weak-coup- 
ling constant  scheme. 
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